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Abstract 

We solve the approximate renormalisation group found by Qiu Niu and Franco 
Nori Q for a quasiperiodic tight-binding hamiltonian on the Fibonacci chain. This 
enables us to characterize analytically the spectral properties of this model. 
PACS numbers: 05.45+b 61.44+p 71.25c 
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In the past few years, many models |T|]have shown that a particle moving on a infinite 
chain and subjected to a quasiperiodic (QC) or incommensurate (IC) modulation can ex- 
hibit critical localization properties. Depending on the QC (IC) modulation strength, in 
contrast to the case of a disordered modulation, the particle wavefunction is not necessarily 
strongly localized (Insulating); it can also be extended (Metallic) as in a translationally in- 
variant system or in a critical localization regime at the Metal-Insulator transition. Typical 
models where these peculiar localization properties occur are Harper like models (IC) or 
tight-binding hamiltonians associated with a quasiperiodic sequence such as the Fibonacci 
sequence (QC). In general, these models have two important parameters: an irrational uj 
which is responsible for the absence of periodicity, and a parameter K which determines the 
strength of the (QC) or (IC) modulation. The usual procedure to characterize the localiza- 
tion properties of these hamiltonians H w K is to study the scaling of the spectral properties 
of a sequence of periodic hamiltonians H^^k which converge to H Ut K as n goes to infinity. 
If the widths of the q n individual bands that compose the spectrum of H Un K , decrease ex- 
ponentially with the period (q n ) of H UJn x then, H Wj k is i n a insulating regime and has a 
Pure Point spectrum. In the case when the widths decrease only inversely proportionally to 
q n , then H Wt K is in the metallic regime and has a Absolutely Continuous (AC) spectrum. In 
contrast with these two cases, the critical localization regime of Hu,k is expected to have a 
Singular Continuous (SC) spectrum with multifractal (MF) properties in the spectral mea- 
sure. In fact, in this latter case, the scaling found by many numerical simulations [0J is 
that the bandwidths decrease like ~ q^ 1 ^" where the exponent a(< 1) varies from band 
to band so that there are typically ~ q^ 9 ^ bands with the same exponent. Due to these 
peculiarities, the critical regime has been extensively studied by both numerical simulations 
and analytic techniques However, as yet, a analytical quantitative determination of the 
exponents ,a and g(a), has not been achievable for even one irrational uj. 

The purpose of this paper is to show a QC tight-binding model for which we are able to 
analytically characterize all the spectral properties. Our work starts from the approximate 
renormalisation group (RG) found by Niu and Nori for a QC model on the Fibonacci 
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chain. By reformulating and solving the RG of Niu and Nori, we derive constructive and 
transparent recurrence schemes for both the energy levels and the bandwidths. From these 
two schemes we deduce new recurrence relations for the spectral measure, the large time 
average return probability of particle defined in ||, the spectrum Lebesgue measure, the 
MF partition function || and the bandwidth distribution. For most of these relations, a 
natural fixed point solution is a power law, either in size or time. By comparing with the 
fixed point equation of the MF partition function, it appears that the exponents associated 
to these power laws, are related to a subset of the anomalous dimensions which characterize 
the MF properties of the spectral measure. A direct calculation of the function g(a) vs a 
confirms these MF properties. To complete the analysis of the spectral properties, we also 
study the gap properties. We find that there are two types of gaps : transient and stable. 
For the first type, their properties are like those of the bandwidths. In contrast, for the 
stable gaps, the distribution of their width (g) is a stable power law P(g) ~ g~( 1+DF ) where 
Dp is the Hausdorff dimension of the spectrum measure. 

These results complete and correct a previous MF analysis of the spectral measure of this 
model ||; our work also unifies many partial results obtained by other methods, for both 
this model and the Harper model ||. 

We consider the tight-binding hamiltonian H n defined on approximant of period F n of 
the Fibonacci chain by the following equation: 

F " . 

H n = 22 ViC { Ci + tj 5 j + iCjQ + i + tj_i,iCjCj_i (1) 
i=i 

The on site potential Vi is taken to be uniform (Vi = V). In contrast, the hopping amplitude 
t i)i+ i from site i to site i+1 is given by i^+i = t w (l—x{u n i))+t s x{oJ n i)-) where u n = ^f^- tends 
to the golden mean u= v ^~ 1 (F n+ i = F n +F n _ 1 and F n ~ to 11 ). The characteristic function 
x(u;i) takes the value or 1 according to the Fibonacci sequence and correspondingly, the 
bond will take the value t w or t s . For finite n the density of weak bonds (t w ) is u n and 
tends to uj in the quasiperiodic limit. For the strong bonds (t s ) the density is uj\ = -^F^ 
and tends to uj 2 . The periodicity of the hamiltonian H n (V,t s ,t w ) allows us to define Bloch 
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boundary conditions of the form Cj + F n = e Cj. For a fixed k, the energy spectrum of 
H n (V,t s ,t w ), which we define as W n (V,t w ,t s ), consists of F n levels E l n {k) (i = 1, ...,F n and 
E l n < E 1 ^ 1 by convention). Varying k from to 7r allows the association of an energy band 
of width A l n = \E l n (ir) — E l n (0)\ to each of these levels. 

Using a perturbative approach, Niu and Nori have shown that in the strong modu- 
lation regime (t w /t s <C 1) the spectrum W n (0,t w ,t s ) is the union of three sub-spectra 
W n - 2 (V + ,t+,t+), W / n _ 3 (V°,t°,t° ) and W n - 2 (V~ , t~ , t~), which correspond to three sub- 
hamiltonians with periods F n _ 2 , F n _ 3 , F n _ 2 and renormalized parameters V ±,Q , tf'°,t^,' 
respectively. A representation of this perturbative RG, with the explicit value of the renor- 
malized parameters, is schematically given by the following relation (0). 



In principle, the scheme (§) simplifies the problem, since it relates the spectral properties of a 
hamiltonian of period F n to those of three sub-hamiltonians of smaller period. However, it is 
clear that upon I iterations of (0), the difficulty which is initially due to the large period F n , is 
replaced by the problem of an increasing number (3') of different hamiltonians to be solved. 
As we now describe, there is some properties of both the hamiltonian and the RG (0) allow us 
to overrule this difficulty Firstly, it is clear that the spectrum of H n (0, t s , t w ) is independent 
of the sign of t s and t w , thus we have W n (0,t s ,t w ) = W n (0, \t s \, \t w \) |L0]. Secondly, we see 
that the spectrum W n -2{V ± , tf , t^) is just uniformly translated from W n _2(0, tf, t^) by a 
factor V ± . Thirdly, the renormalized parameters have the following property: 



W n {0,t s ,t w ) 




(2) 



tf\=zt 



s ) 



Z = 




w 



(3) 



Combining these properties and using (Q), we deduce the following new renormalisation 



scheme: 
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-t s + zW n - 2 {0,t s ,t w ) 
W n (0,t 8 ,t w ) — > zW n . 3 (0,t s ,t w ) (4) 
+t s + zW n - 2 (0,t s ,t w ) 

We see that if t s is sufficiently strong, the spectrum W / n _ 2 (0, t s , t w ), which is contracted by a 
factor z and centered around ±t a , is not mixed with the spectrum W n _ 3 (0, t s , t w ), which is 
contracted by z. More precisely, if we call A n the distance between the lowest and highest 
energy levels of H n , then, this non-overlapping condition becomes (zA n _ 2 + zA n s) < 2t s . 
Under this condition, relation (|j) gives the following recurrence scheme between the energy 
levels Ei(0,t s ,t w ), ^_ 2 (0, t s , t w ) and ^_ 3 (0, t„ t w ): 

E l n =-t A + zE i n _ 2 (i = l,F n _ 2 ) 

E i+ Fn _ 2 =zEU (< = l,F B _s) (5) 

F^ +Fn_1 =t A + zEi_ 2 (i = l,F„_ 2 ) 
Similarly, the associated recurrence for the bandwidths is given by: 

An = 

A j+F n _ 2 = (6) 

We now give the quantitative consequences of the last three relations (§)(§]) (@) upon 
the spectral properties. From the recurrence scheme (|5|) we see that we can assign a set of 
indices {+}(+t s ), {0} or { — }(— t s ) to each individual energy level, according to the path of 
bifurcations of that level. Therefore a typical level E l n has n + , Uq and n_ indices (+, 0, — ); 
with the constraint 2(n + + nJ) + 3no = n(±l). As examples, the lowest level is indexed by 

{ ...}, (n_ = [n/3],n + = n = 0) and the highest is indexed by {+ + +...} [0 |Tl| . From 

this indexation and relation @, we see that the band associated with a level E % n (n +) no, n_) 
has a width A l n (p, q) pd z p z q with p = (n^+ri- ) and q = riQ. Consequently, the number of bands 
of width A n (p, q) is given by N n (p, q) = 2 P {^ +q ^j- These last two results allow us to calculate 
the exponents a and g(a) defined in the introduction (A^ = F~ l l ai and N n (p, q) = F^ a ^). 
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As shown in relation (0), in the quasiperiodic limit (n — > oo), these exponents are function 
of the parameter x = p/n which varies continuously in [0, 1/2] 0. 

a(x) = \nuj/(x\nz/z 2 / 3 + lnz 1 / 3 ) 

<7(a(ar)) = (x In 3x/2 - (1 + x) In (1 + xf' 3 (7) 
+ (1 - 2x) ln(l - 2x) 1 / 3 )/lnw 

From this last relation we can obtain two interesting properties of the spectrum. Firstly, 
we observe that when z = z 2 / 3 {t w /t s = 1/8), the exponent a is independent of x. As a 
consequence the spectrum is a pure fractal with Hausdorff dimension Dp = a = \nuj 3 /liiz. 
In contrast, when for example z > z 2 ^ 3 , we see that the exponent a varies between a m i n = 
lnuj 3 /\nz and a max = In uj 2 / In z. As these two exponents correspond to the bandwidth 
associated with the central and edge levels respectively, this property allows us to compare 
their values with the exact analytical result of Kohmoto [[7] H.From this we observe that our 
shrinking factors z and z are just the dominant terms in a (t w /t s ) series expansion of the two 
more exact values z ex = 2/(^/(7-1)2-4 + 7-1) and z ex = l/(^/l + 4(7 + l) 2 + 2(1 + 1)) 
where I = U¥ - r^) 2 and J = 3 + V25 + 167. 

4 v t s t w ' 

So far, we have only analyzed the properties of individual bandwidths and levels. However 
it is also possible and very instructive to study integrated quantities. We start with the 
spectral measure and a physical quantity closely related to it. The spectral measure dix n (E) 
of the hamiltonian H n is dfx n (E) = p n (E)dE where p n (E) = y- Yn=i ^{E — E 1 ^) is the density 
of states. From this definition and relation (||D we deduce the following recurrence || : 

dfi n (E) = LU 2 n dp n - 2 { s ~) + ^d/i n _ 3 (— ) + L0 2 dn n -2{ s ~) (8) 

z ' z z 

As an example of the application of relation ([D, we calculate the large time average return 
probability defined by p n (t) = | /+~ e~ iEt d[i n {E)\ 2 = 2njl(t)fi*(t) §. Using © we firstly 
deduce that jl n {t) = 2u 2 cost s t/2 n _ 2 (^t) + uj 3 fl n _ 3 (zt) . Now, in the large time limit, we have 
on average (cost s t) ~ and (cos 2 t s t) ~ 1/2, and from this we immediately get: 

p n (t) = 2^ n p n _ 2 {zt) + u%p n - z (zb) (9) 
6 



In the limit (n — > oo), we see that a fixed point solution of relation (0) is ~ t~ 7 where 
the exponent 7 is determined by 2oj 4l z~ 1 + a; 6 z _7 = 1. As we now show, this exponent 7 
is one of the anomalous dimensions D q that characterize the MF properties of the spectral 
measure. In our case, these non trivial dimensions D q are defined by the requirement that 
the partition function [f| T n (q, r = (q — l)D q ) = F~ g Y^i=i(A t n )~ T De stationary in the limit 
n — > 00. Using relation (J6|) we get the following recurrence for the T n (g, r) ||: 

T n (g, r) = 2^r„_ 2 (g, r) + ^r n _ 3 (g, r) (10) 

The stationary constraint then gives a self-consistent equation for the D q 

2uJ 2 qz (.l- q )D q+ ^ q -{l- q) D q = l (n) 

From this last relation, we immediately see that the exponent 7 previously defined is in fact 
equal to D2. A second consequence of (|TID is that the Hausdorff dimension Dp = Dq is the 
solution of 2z Df + z Dp = 1. To see further the use of relation (|TT| ) and the role of the D q 
we calculate two other quantities of interest: the Lebesgue measure B n = J2i=i ^ni an d the 
number of bands of width between A and A + dA, dN n (A) = Yh=i 5{A — A l n )dA. Using 
relation (^) we can easily deduce a recurrence relation for each of these quantities: 

B n = 2zB n _ 2 + zB n ^ (12) 

dN n (A) = 2dN n ^) + dN n _ 3 (j) (13) 

The first of these equations was partially guessed in for a model on the Fibonacci chain; a 
very similar relation was also derived for the case of the Harper model in 0. From equation 
( 0), we can show that the large n behavior of B n is B n ~ B F~ S where the exponent 
6 is related to the anomalous dimensions by D_$ = y+5- Now, considering ( [L3|), we see 
that in the limit n — > 00, a possible invariant form is given by dN*(A) = A-^+^dA with 
(3 = Dp- In view of relation (0), this simple invariant solution is quite surprising and indeed 
its sense is not very clear. In particular, if instead of dN n (A) we look at the distribution of 



bandwidths W n {A) = d p n ^ , then the corresponding recurrence relation allows the two 
possible invariant distributions W*(A) = 6(A) and VT*(A) = A -1 , whose form do not 



coincide with W*(A) = = A^ 1+D ^ @ 



To complete the study of the spectral properties, we also look at the the statistical 
properties of the gaps. Roughly speaking a gapwidth is the distance between two levels. In 
consequence, we might expect their distribution to be quite similar to that of the bandwidths. 
However, there is an important difference; the number of gaps of a chain F n is F n — 1, thus, 
2(.F n _2 — 1) + F n s — 1 = F n — 3 < F n — 1! Looking at figure [I] the last inequality means that 
if we take only gaps coming from zW n -2(±t s , t s , t w ) and zW n -z(Q, t s , t w ), we miss two gaps 
which are in fact precisely the biggest. Taking this into account, we see that the number 
of gaps of width between g and g + dg, dN n (g) = Y^i=x${9 — 9n)dg, obeys the following 
recurrence (n > 3): 

dN n (g) = 2dN n . 2 (^) + rfiV n _ 3 (|) + 26{g - g° n )dg (14) 

where the last term re-introduces the two largest gaps of width g>° at each iteration. A further 
refining of our description requires two other important remarks. The first is that the initial 
conditions are dN (g) = 0, dN\(g) = and dN 2 (g) = S(g — go)dg. The second is that in the 
limit n — > oo the width g® tends to a fixed value g* = t^ — \{zA* + z A*) = t& ~wE^ where 
A* is the width of the spectrum in this limit. This last remark allows to replace the previous 
recurrence for dN n (g) by the following effective equation which describes the infinite size 
behavior more effectively: 

dN n (g) ~ 2dN n ^) + rfiV n _3(|) + 25(g - g*)dg (15) 

As shown in table | the iteration of relation (|T5|), with the previous initial conditions, pro- 
duces two kinds of gaps. The first kind corresponds to what we call the transient gaps. 
For a chain F n , these gaps have widths of the form g = z p z q go with (2p + 3q = n — 2) and 
N n (p, q) = 2 p (p +l? ) (last two columns of table [I]). These transient gaps are created by iter- 
ation of the initial condition dN 2 {g) = S(g — go)dg and their effective recurrence does not 
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contain the last term of flip]). For these reasons their distribution is strictly similar to that of 
bandwidths and in particular the sums of their widths decreases like the spectrum Lebesgue 
measure B n . In contrast, the second kind of gaps are those created by the presence of the 
last term in relation (first two columns of table | ). As can be seen in table |, if a gap 
of this kind opens for, say a chain F p , it persists for longer chains; it is stable. For a such a 
chain F n , these stable gaps have widths of the form g = z p z q g* with N n (p } q) = 2 P {^ +q ^j but 
now 2p + 3q takes all values between and n — 3. Due to this difference the distribution 
of the stable gaps differs from that of the transient gaps in the following way: in a similar 
manner to ([13]) , the absence of the last term yields an invariant solution to (|15|) of the 
form dN*(g) = g~( 1+£>F ^dg. However, similarely to bandwidths |T!|, for the transient gaps, 
dN n (g) tends to dN*(g). In contrast when the last term is present, that is for the stable 
gaps, then, the function dN n (g) really tends to dN*(g) = g~( 1+Dp ^dg over the whole interval 
[0,(7*]. Due to this propertyin that case, we can also define a distribution which is of the 
form P*(g) = dN d g 9 " 1 — g~( 1+£>F \ An additional property of the stable gaps concerns the 
value (G n ) of the sum of their widths. From ( |T5| ) we see that G n obeys a recurrence relation 
G n = 1zG n -2 + zG n -3 + 2g*, and from this we can deduce that in contrast to B n , G n does 
not decrease with F n but tends to a value G* which is exactly equal to the spectrum width 
A*. 

In conclusion, we have described as completely as possible the statistical properties of 
the energy spectrum of a tight-binding hamiltonian on the Fibonacci chain. We have com- 
pared several of the new predictions with numerical computations with satisfactory results. 
As the text has made clear, our results, qualitatively and quantitatively complete and unify 
previous works on similar models. 
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TABLES 

TABLE I. Gapwidths obtained for the six first iterations of relation (15). Stable gaps are 
written in the first column. Transient gaps are in the third column. 
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FIGURES 

FIG. 1. Spectrum of the approximant hamiltonian (H n ,F n ), n < 8 deduced from relation (5) 
for [t s = l,t w = 0.5). The three initial conditions are, (Ho, Fq = l,ti = t s ); (H±, F± = 1,U = t w ) and 
(H2,F2 = 2,t2i =t s ,t2i+i =t w ). z and z are the two shrinking factors, go and g* , are the initial 
transient gap and the maximum stable gap respectively. 
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